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We compute the complete topological expansion of the formal hermitian two- 



' matrix model. For this, we refine the previously formulated diagrammatic rules 

for computing the i expansion of the nonmixed correlation functions and give a 
new formulation of the spectral curve. We extend these rules obtaining a closed 



formula for correlation functions in all orders of topological expansion. We then 
integrate it to obtain the free energy in terms of residues on the associated 
Riemann surface. 
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> : 1 Introduction 



Matrix models is a fascinating topic unifying many otherwise seemingly unrelated 
disciplines, including integrable systems, conformal theories, topological expansion, 
etc jSHlEl- However, only recently the understanding came that geometry itself plays 
a crucial role in constructing perturbative solutions to matrix models. Dijkgraaf and 
Vafa ^H] studies in the one-matrix model (1MM), together with subsequent observa- 
tions [2B1IHIHDIH2] based on the idea of [SOj that topological hierarchies are encoded 
already in the planar limit of the solutions to the Hermitian one-matrix model with 
multiple-connected support of eigenvalues, led to understanding of a crucial role of the 
corresponding spectral curve in constructing explicit solutions. 

The first attempts of computing the sub leading terms in the 1/N 2 expansion, where 
based on loop equations [2113] , and allowed the authors of [3] to find the first few terms 
of the expansion of the free energy, only in the 1-cut case of the 1-hermitean matrix 
model. 
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In the begining of the 2000 's, a geometrical approach, supported by using the 
master loop equation, resulted in constructing solutions in the first subleading order 
of the 1MM in the mutlticut case HSHnjCSHESl • 

Then, a Feynmann-like diagrammatic technique was invented in ^H]; which allowed 
to reformulate the loop equations themselves in a proper geometric way, and it became 
possible to compute all correlation functions of the 1-MM to all orders in the 1/N 2 
expansion [T^]. However, the diagrammatic technique of JH] could not be applied 
directly to the expansion of the free energy. The topological expansion of the free 
energy to all orders was found with a refined diagrammatic technique in [TT], where 
the key ingredient was the homogenity property of the free energy to all orders. 

In parallel, the hermitian two-matrix model (2MM) solutions have been obtained, 
at early stages, in the planar limit of the 1/iV-expansion |Mll21j . It was however almost 
immediately observed jSlCHllS] that the 2MM solution in the planar limit enjoys the 
same geometrical properties as the 1MM solution; only the spectral curve becomes 
an arbitrary algebraic curve, not just an hyperellitic curve arising in the 1MM case. 
The subsequent progress was however hindered by that the corresponding master loop 
equation in the 2MM case cannot immediately be expressed in terms of correlation 
functions alone. Nevertheless, using the geometrical properties of this equation, the 
solution in the first subleading order has been constructed [1211201 on the base of 
knowledge of Bergmann tau function on Hurwitz spaces. However, the general belief 
arose that the 2MM case should not be very much different from the 1MM case, that 
is, we expect to find all ingredients of the 1MM solution in the 2MM case. Next step 
towards constructing the topological expansion of the 2MM was performed in [21], 
where the first variant of the diagrammatic technique for the correlation functions in 
the 2MM case was constructed. In the present paper we improve the technique of [23 
(actually effectively simplifying it) to accommodate the action of the loop insertion 
operator. In fact, we demonstrate that the diagrammatic technique for the 2MM case 
is even closer to the one in the 1MM case that was before: in particular, we need only 
three-valent vertices, and the additional operator H we need to obtain the free energy 
turns out to be of the same origin as the one in the 1MM. 

The paper is organized as follows. In Sec. 2, we collect all the definitions, algebraic- 
geometrical notation, and facts about loop equations and filling fractions we need in 
what follows. In Sec. 3, we provide a new formula for the spectral curve, which results 
in the new diagrammatic rules formulated in Sec. 4. In the same section, we express 
the action of the loop insertion operator in terms of our diagrammatic technique, which 
makes the construction closed as regarding the nonmixed multipoint correlation func- 
tions. We introduce the "integration" operator H in Sec. 5 and, using this operator, 
present the diagrammatic technique that enables us to construct the complete topologi- 
cal expansion for the 2MM to all orders except the subleading order. But for this latter, 
the answer has been found in [20J, so we eventually formulate a complete procedure 
for constructing the topological expansion in the 2MM. In Appendix A, we extend our 
technique to calculating the first mixed correlation function; higher mixed correlation 
functions need further refining of this technique, which is beyond the scope of this 
publication. In Appendix B, we prove the symmetricity of the free-energy expression 
w.r.t. interchanging x and y variables, which is crucial for the proper integration of 
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the expression for the first mixed correlation function. 

2 Definitions and algebraic-geometrical notation 

2.1 Definition of the model 

We study the formal two-matrix model [27] and compute the free energy T of this 
model in the asymptotic -^-expansion. The partition function Z is the formal matrix 
integral 

Z:= f dM 1 dM 2 e-^KMMO+W)-^) = (2 _ x) 

J Hn x Hm 

where M\ and M2 are two N xN Hermitian matrices, dM\ and dM.2 are the products of 
Lebesgue measures of the real components of M\ and M 2 , h = ^ is a formal expansion 
parameter and Vi and V2 are two polynomial potentials of respective degrees d\ + 1 
and <i 2 + 1 

^) = £^ fc , %(y) = ^|y*. (2-2) 
fc=i fc=i 

Formal integral means that it is computed order by order in powers of the t^s (see 
section l2~3*l or |2*2*]). We consider polynomial potentials here only for simplicity of 
notations, but it is clear that the whole method presented in this paper extends to 
"rational potentials" (i.e. such that V{ and V 2 ' are rational functions), and we expect 
it to extend to the whole semiclassical setting [6 including hard edges as well. 

One can expand the free energy as well as all the correlation functions of this model 
in a h series assuming h to the order of the reciprocal matrix size; this procedure is 
customarily called the topological expansion pertaining to the fat-graph representation 
for formal integrals ( jHSHHllIl] ) . 

The topological expansion of the Feynman diagrams series reads, in terms of the 
free energy, 

00 

T = F(N, T, f ls t 2 , . . . , t dl+1 , h, h, ■ ■ • , U 2+1 ) = h 2h - 2 ^ {h) (T, ...). (2-3) 

h=o 

We find a general formula for J 7 ^ 1 ' for any positive integer h. Actually, we address 
the problem for h > 2; the solutions for h = ( [HIHIl]) and h = 1 ( 20 J are already 
known. 

2.2 Notations 
2.2.1 Variable sets 

We consider functions of many variables x\, X2, x%, . . ., or of a subset of those variables. 
For this, we introduce the following notation: 
Let K be a. set of k integers: 

K = (i 1 ,i 2 ,...,i k ). (2-4) 
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Let k = \K\ denote the length (or cardinality) of K. For any j < \K\, let Kj denote 
the set of all j — upples (i.e., subsets of length j) contained in K: 

K 3 :={JdK , \J\=j}. (2-5) 

We define the following k— upple of complex numbers: 

k 

Xk := (x h ,x i2 , . . . ,x ik ) , dx K :=Y\_dxj. (2-6) 

3=1 

2.2.2 Resolvents 

For given integers k and I, we define the resolvents: 

. 9 <9 d d d dT 

= —n 



conn 



dV 2 {y\ L \) • • • 9F 2 ( yi ) dV l {x m ) dV l {x\ K ^ l ) ' ' ' dV^m) 

+ WW + Vl,0^|L|,l 



»i Vi 

(2-7) 



with the formal loop insertion operators 



3 9 a 9 sr J 9 to &\ 



9Vi{x) j^x^dtj dV 2 (x) j^y j+1 dtj 

We also introduce the polynomials in y: 

\ r — 1 / conn 

where the subscript conn denotes the connected component. 
For convenience, we renormalize the two-point functions: 

w\k\,\l\(xk, Yl) = w\ Km (x K , y L ) + _ + _ (2-10) 

and their polynomial correspondent: 

u k (x,y;x K ) ■■= u k (x,y;x K ) - 5 k:0 (V 2 (y) - x). (2-11) 
We consider the ^-expansions of the above quantities: 



w\k\,\l\(xk, Yl) = J2 ^ w \k\,\l\(*k> yi)> ( 2 " 12 ) 

h=0 



u k (x,y;x K ) = y^h 2h uf\x,y : x K )- (2" 13 ) 

h=0 

We also need the functions closely related to the algebraic structure of the problem: 

Y(x):=Vl(x)-w 1 (x) (2-14) 

and 

f|r ,, ): . t / t ,!M I gM)\ . (2 . 15) 

\ 1 o & l conn 

The latter function together with all its terms of ^-expansion, is a polynomial of 
degree d\ — 1 in x and d 2 — 1 in y. 

2.2.3 The master loop equation 

Among many different ways to solve matrix models [33J, addressing the formal 2- matrix 
model problem, we choose the so-called loop equations [SUE] encoding the matrix 
integral in ()2-lj) to be invariant under special changes of variables. They correspond 
to the generalization of the Virasoro constraints of the one-matrix model, i.e. the 
W-algebra. 

Considering a particular change of variables, we come to the master loop equa- 
tion Pl|: 



E(x, y) = (y- Y(x))u (x, y) + f?u x {x, y; x) 



(2-16) 



where 

E(x,y) = (V{(x) - y){V 2 \y) - x) - P(x,y) + T (2-17) 

is a polynomial in x and y defining the spectral curve. 

Considering the 't Hooft expansion in orders of h 2 , for any h > 1, we have: 

E {h \x,y) = (y- Y(x))vg»(x,y) + w[%)u^(x,y) 

hl d 

m=l ^ ' 

(2 - 18) 

where E^ h \x^y) is the /i'th term in the ^-expansion of the spectral curve E(x,y). 
2.2.4 Algebraic geometry notations 

In this section, we recall the algebraic-geometrical pattern of our problem (see more 
details in [211213123] ) • 

To leading order, the master loop equation (|2-18|) reduces to an algebraic equation 

E {0) (x,Y(x)) = (2-19) 
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with E^°\x, y) being a polynomial of degree &\ + 1 in x and d 2 + 1 in y. 

We parameterize the algebraic curve y) = implied by (15131) by a running 

point p on the corresponding compact Riemann surface £. We therefore define two 
analytical meromorphic functions x(p) and y(p) on £ such that: 

E {0) {x,y) =0 ^3\p e £ x = x{p) , y = y(p). (2-20) 

The functions x and y are not bijective. Indeed, since E^°\x, y) has a degree rf 2 + 1 in 
y, it admits d 2 + 1 solutions for a given x; that is we have d 2 + 1 points p on £ such 
that x(p) = x. Thus, the Riemann surface is made of d 2 + 1 x-sheets, or respectively, 
of g?i + 1 y-sheets. We then denote 

x(p) = x 4=> p = p^\x) for j = 0, . . . , c?2, (2-21) 

y(p) = y p = pV\x) for j = 0, . . . , d%. (2-22) 

Among the different x-sheets (resp. y-sheet), there exists only one where 
yip) ~x( P )^oo V{(x(p))-^ + 0(x 2 (p)) (resp. x(p) ~ y(pH oo V 2 \y(p)) - + 0(y 2 (p))). 
We call it the physical sheet, and it bears the superscript 0. 

Genus and cycles. The curve £ is a compact Riemann surface of finite genus 
g < d\d 2 — l 1 . We do not require g to be equal to d x d 2 — 1, that is, we allow double 
points on the corresponding Riemann surface. We choose 2g canonical cycles as Ai, 
Bi, i — 1, . . . , g, such that: 

Ai n Aj = , Bi n Bj = , AiH Bj = 5ij. (2-23) 

Branch points. The x-branch points fi a , oi = 1, . . . , d 2 + 1 + 2g, are the zeroes 
of the differential dx, respectively, the y-branch points up, (3 — 1, . . . , d\ + 1 + 2g, are 
the zeroes of dy. We assume here that all branch points are simple and distinct. Note 
also that Ey°\x(p),y(p)) vanishes (simple zeroes) at the branch points (it vanishes in 
other points as well). 

A branch point is a point where two sheets of the Riemann surface meet. If p is a 
point on the Riemann surface near a branch point, there is another point p l , which we 
note p, near the same branch point. In other words: 

Va 3\p 7^ p such that x(p) = x(p) and y{j>) ^ p ^^ a y{p) (2-24) 

Let us emphasize that the point p and the corresponding sheet, depend on which 
branch point we are considering. 

Bergmann kernel. On the Riemann surface £, we have a unique Abelian bilinear 
differential B(p, q), with one double pole at p = q such that 

B{p,q) - d ^P) dx (l) + finite and v , I B (p,q) = 0. (2-25) 

P^q Kg) - x(W 2 JpeAt 

1 This genus g must not be confused with the genus indicating the order of the topological expansion 
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It is symmetric, 

B(p,q) = B(q,p), (2-26) 

it can be expressed in terms of theta- functions [2*6ll25| . and depends only on the complex 
structure of £. 

Abelian differential of the third kind. 

On the Riemann surface S , there exists a unique Abelian differential of the third 
kind dS qjr (p), with two simple poles at p = q and at p — r, such that 

Res dS qr (p) = 1 = — Res dS qr (p) and Vi (p dS qr (p) = 0. (2-27) 

p^q ' p^r ' _/ A 

Notice that the Abelian differential of the third kind and the Bergmann kernel are 
linked by 

dS q , r (p)= f 9 B(£,p) and B(p,q) = d q (dS q>r (p)) . (2-28) 

J£=r 

where the contour of integration is a line from r to q, which does not cross any A or B 
cycle. 

Given a branch point p a , and a point q in the vicinity of /i a , we introduce the 
following notation: 

dE qrq (p)= Fb(£, P ), (2-29) 

Jq 

where now, the integration path is chosen as the shortest path between q and q, i.e. a 
path which lies in a small vicinity of p a . That definition differs from the one above. If 
the branchpoint p a is surrounded by contour Ai, we have: 

dE qrq (p) = dS qrq (p) + <f B(p, (2-30) 

The main property of that dE q ^ q (p), is that it vanishes at q — q, i.e. it vanishes at the 
branch point p a . 

Correlation functions on the Riemann surface 

Given the algebraic curve, we see that we can redefine the correlation functions 
fl2-12|) and (J2-13|) more precisely. Indeed, they are defined only as formal series as their 
arguments and tend to infinity. The loop equations show that these formal series 
are in fact algebraic functions and hence have a finite radius of convergency, with cuts 
beyond the radius. As functions of the x and y variables, they are multivalued. 

If instead of writing them as functions of x or y, we write them as functions on the 
Riemann surface, they become monovalued. This is the reason why we prefer to in- 
troduce the following notation for the correlation functions, as meromophic differential 
forms on the Riemann surface: 

W\ K \,\l\(pk,cu) ■=w\ K \ t \ L \{x{T£t K ) 1 y{(\L))dx{ii K )dy{(\L), (2-31) 

and 

U\k\(p,v;Pk) :=u\ K \(x(jp),y\x(jpK))dx{p K ), (2-32) 
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where p^s and g/s are points on the surface £ whose images x(pi) or y(qi) are complex 
numbers, y is a complex number. 

Remark 2.1 The interpretation of those correlation functions as generating series for the 
moments < Yli Tr M^' Y[j Tr M 2 J > c , corresponds to the situation where all p^s are in the 
x-physical sheet, in a vicinity of oo x , and all q^s are in the y-physical sheet in a vicinity of 

OOy. 



2 point function 

It is well known ( see for instance P^ l3lHl2^ l5l l2^ll9| l2T| ) that, to leading order, the 
2-point function W^o is the Bergmann kernel: 

W^(p,q) = B(p,q) (2-33) 
The non-renormalized 2 point function: 

W 2Ap , q) = W 2 , oM - J^L (2-34) 

is finite at p — q. 



2.3 Loop equations and fixed filling fractions 

We have showed that to leading order, the 1-point function Y^(x) obeys an algebraic 
equation (|2-19|) : 

E^\x,Y{x)) = 0, (2-35) 

where 

E^(x,y) = (V((x) - y){Vi{y) - x) - pW(x,y)+T (2-36) 

But so far we have not discussed how to determine the polynomial P^ (x, y) . 

As this was extensively discussed in the literature, we only briefly summarize it 
below. 

We need did?, — 1 equations to fix the d\di — 1 unknown coefficients of p(°\ Those 
additional d\di — 1 equations do not come from any loop equation, so we need an 
independent hypothesis related to the precise definition of our matrix model. In a 
sense, loop equations express the invariance of the integral under reparameterizations, 
independently on the integration paths. Additional equations are those that depend 
on the choice of integration path. 

For arbitrary integration paths, the h expansion may not exist. The choice of 
integration method we use in this paper, corresponds to the so-called formal matrix 
model, with fixed filling fractions. 

We consider here the problem of a formal matrix model, i.e. the one represented by 
the formal power series expansion of a matrix integral, where the non-quadratic terms 
in the potentials V± and Vi are treated as perturbations near quadratic potentials. Such 
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a perturbative expansion can be performed only near local extrema of Vi(x)+V 2 (y)—xy, 
i.e. near the points (£j, rji), i — 1, . . . , did 2 , such that 

n'(&)=^ and V 2 \ Vl )=^ (2-37) 

which has in general did 2 solutions. Therefore, perturbative expansion can be per- 
formed near matrices of the form: 

ni U2 n did 2 



Mi = diag fi, &»•••» 6» ••• » • -T&mJ, ( 2 - 38 ) 

M 2 = diag (f/77~^~^i , ^27^~^2, • • • , ^did 2 , • - ,%d 2 v ) , ( 2 - 39 ) 

such that Xlfii 2 n « = The perturbative integral is computed by writing M\ = 
Mi + 5 Mi and M 2 = M 2 + 8M 2 , by expanding higher order terms (cubic and higher) in 
5 Mi and 5M 2 , and by treating quadratic terms in 5 Mi and 5M 2 as gaussian integrals. 

In other words, a formal integral can be computed as soon as we have chosen a 
vacuum around which to make a perturbative expansion. The choice of a vacuum is 
equivalent to choosing a partition of iV into d\d 2 parts: 

d\d2 

J2ni = N. (2-40) 

i=i 

It is easy to see that if we truncate the perturbative expansion to any finite order, 
we have: 

— I ( tr \ dx = — I tr l -=dx = -m, (2-41) 

2i7cf Ci \ x-Mil 2i7cf Ci x-Mi K J 

where Cj is a small direct circle around the point £j in the complex plane. And thus: 

Wifi(x)dx = —riih := — q. (2-42) 



1 

27tt 



If we do not truncate to a finite order, all functions become algebraic, with cuts, and 
the contours C« are then enhanced to finite contours around the cuts, and up to a 
redefinition of the contours, the filling fractions are the .4-cycle integrals. We define: 

iii l vdx= -^i xdy=e ' (2 - 43) 

We call ej's the filling fractions, and they are given new parameters (moduli) of the 
model. In what follows, we consider them to be independent of the potential or on any 
other parameter. Let us notice that the such defined q's are not necessarily of the form 
Ui/N with Hi a positive integer, but they can be arbitrary complex numbers, provided 
that: 

e i = t = Nh (2-44) 
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In particular, since all correlation functions Wk,o(xi, ■ ■ ■ , Xk) are obtained by deriva- 
tion of Wifi with respect to the potential Vi, we have: 

7p- / w k>0 (xi, . . . ,x k )dxi = 0. (2-45) 

Eq. (j2-43|) together with the large x and y behaviors, suffices for determining com- 
pletely all the coefficients of the polynomial E^°'(x,y), and thus the leading large N 
resolvents wi t o(x) and Wo,i(y)- 

Remark 2.2 It is easy to see that if we truncate the perturbative gaussian integral to 
any finite order, the result is a polynomial in 1/N 2 , and thus, formal matrix models always 
have a 1/N 2 expansion. This comes from the fact that, to any finite order, the perturbative 
gaussian integral is the generating function for counting discrete surfaces made of a finite 
bounded number of polygons, thus having a maximal Euler characteristic (power of N). 

Remark 2.3 We emphasize again the formality of our model: in principle, the ^4-cycles do 
not necessarily lie in the physical sheet, so we would need additional assumptions to segregate 
physically meaningful models. Recall that the ej's can be arbitrary complex numbers. 

Remark 2.4 We have d\d2 non vanishing contour integrals, but they are not all indepen- 
dent. Their sum can be deformed as a contour around oo. Thus, we only have d\d% — 1 
independent non vanishing contour integrals, and thus the maximal genus of our curve is 
d\di — 1. The number of independent filling fractions is d\d<i — 1, due to the constraint on 
their sum. 



Double points The Riemann surface 8 may have double points. 
A point p G £ is a double point iff: 

3q ^ p, x{p) = x(q) and y{p) = y(q). (2-46) 

Double points are such that both derivatives E x and E y vanish, and thus the curve has 
two tangents at such points. The differential forms dx and dy do not vanish at double 
points. 

We need an extra asumption to deal with double points. 

We may consider that in the framework of the perturbative gaussian matrix integral 
seen above, double points correspond to vanishing filling fractions e$ = 0. Thus: 

jw lfl {x)dx = Q (2-47) 

where C is a contour which encircles the double point. In other words, we make the 
assumption that the resolvent has no residue at double points, to all orders in the 
h expansion. We will see below that this vanishing residue condition is sufficient to 
ensure that correlation functions never have poles at double points. 

Another reason to make that assumption is based on the Krichever formula |3Tj] , 
which involves summation only over branch-points and not over double points. 
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3 A new formula for the spectral curve 



The 't Hooft expansion of the master loop equation (|2-18p . links correlation functions 
of genus h to correlation functions of lower genus. The knowledge of the LHS E(x,y) 
would then give a way of deriving the whole topological expansion of the W\k\, o- 
From the definition of the surface £ , one derives the leading term as: 



<ii 



E^(x(p),y(q)) = -t dl+1 xH(x(p)-x(q®(y))) 
= -t d2+1 xH(y(q)-y(pV(x))) 



(3-1) 



i=0 



Actually, not only the leading term, but the whole function admits such a product 
structure: 

Theorem 3.1 The functions E(x,y) and Uo(p,y) can be written: 



d-2 



E(x(p),y) = -t d2+1 " ( l[(y - V{{x(p)) + HTr 



1 



vi=0 



x(p i ) - Mi 



and 



U (p, y) = " ( l[(y - V((x(p)) + hTr 



x(p { )- Mi 



(3-2) 



(3-3) 



where the quotes " (.) " mean that when we expand into cumulants, each time we find 
a 2-point function W2$, we replace it by w 2 ,o = ^2,0 + ( Xl l X2 )2 as ^ n REUB) - 

In other words, E(x,y) is the (c?2 + 1)— point correlation function, with all points 
corresponding to the same x in all the di + 1 sheets, and Uo(p,y) is the rf 2 -point 
correlation function, with all points corresponding to the same x(p) in all sheets but 
the same as p. 
proof: 

We prove this theorem by showing that the two sides of the equalities ()3-2|) and 
(13- 3j) are defined by the same recursion relations. 

We let E(x(p),y) and U (p,y) be defined by the respective RHS of (|3-2|) and 
Consider their topological expansions: 



E(x,y) = J2^ h E {h) (x,y) and U (p, y) = £ H 2h U^(x(p), y) 



(3-4) 



h=0 



h=0 



By expanding in cumulants (connected parts), one can observe that: 
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Et h \x(p),y) = (y - y( P ))Ut\ P ,y) + W^(p)U^(p,y) 

+ £ w^(p)ut m) (P,y) + w^W)^^ 



(3-5) 



m=l 



which has exactly the same form as (J2-18|) . 

Obviously we have E^ \x,y) = E^(x,y) and U^°\p,y) = U^(p,y). 

Let us now show that, given E^°>(x,y), (j2-18|) admits unique solution. For this, 
assume that we know E^ h ~ l \x,y) and UQ m \p,y) for all m < h — 1 and prove that 
(j2-18|) allows computing E^ h \x,y) and Ul^\p,y). 

Consider (j2-18|) for p = q l and y = y(p) for any i = . . . d 2 . It reads 

E^xipUtf)) = £ w[f{f)U { t m \p\y{p 1 )) + Q^r-rft^ip'Mp 1 )^ 1 )- (3-6) 

m=l \P / 

By the recursion hypothesis, one knows the RHS totally, so one knows the LHS. This 
gives the value of E^ h \x, y) for d 2 + 1 values of y. Because E^ h \x, y) is a polynomial in 
y of degree d 2 , an interpolation formula gives its value for any y. It remains to compute 
Uq (p,y) which is straightforward due to (j2-18|) . 

Thus, because E^ h \x,y) and E^(x,y) obey the same equations with the same 
initial condition, we conclude that E^ h \x,y) = E( h \x,y) for any h. 

□ 



4 Diagrammatic rules for the correlation functions: 
new insight. 

Using loop equations, two of the authors have derived two sets of diagrammatic rules 
allowing to compute any WW^o as residues on the Riemann surface [21]. One of these 
theories involved only cubic vertices but presented the disadvantage of using explicitly 
the auxiliary functions U\ K \, whereas the second one involved only the W\ K \pS but 
expressed it through vertices of valence up to d 2 . 

In this section, we present new diagrammatic rules composed of cubic vertices and 
involving only W'a whose arguments are in the vicinity of the branch points. They are 
much simpler than all the preceding ones. 

From ()H-2j) and the definition ()2-17)) of E(x,y), we obtain the following equation: 

-t d2+1 "(]^Uy-V{(x(p)) + hTr^^))- = (V?(y) - x(p))(V{(x(p)) - y) 

-P(x(p),y) + 1. 

(4-1) 

Recall that P(x,y) is a polynomial both in y (of degree d 2 — 1) and in x. 
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The large y expansion of the LHS reads 



i=0 



i=0 j^i x v 



x frTr 



+0(y 



d 2 -2\ 



x(p>)- Mi 



V;(x(p))\)»y 



"^,^2-1 



(4-2) 

and the large y expansion of the RHS yields: 

- t d2+1 y d > +1 - y d \t d2 - i d2+1 V{(x(p))) - t d2+x Q{x{p))y d >- 1 + 0{y d ^ 2 ) : 

where Q(x(p)) is a polynomial in x of degree d\ — 1. 
• Equating the coefficient of y d2 gives 



d-2 



i=0 



l d 2 +l 



which is a polynomial in x(p). Expanded to order h, this means that 

d 2 



X^oVHO for/i>l. 



i=0 



Applying the loop insertion operator to (j4-4j) we get: 



d 2 



E^,o ) (P 4 ^)=^o 



i=0 



dx(p)dx(q) 
x(p) - x(q)f 



which can also be written: 



W 2fi {p\q) + Y,W2fi{p i ,q) = V- 



Equating the coefficient of y d2 1 gives 

d 2 



Crf 2 +1 ^2 + 1 



where 



Q(x) = h < Tr 



V?(s) - V?(Mi) 
x - Mi 
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is a polynomial in x of degree d\ — 1. Using (|4-4|) and (|4-7|) we get: 



d 2 



^ y(p*) 2 + h 2 m 2 , (p\p*) = v[{x(p)f - 2 7 



td 2 -i 



i=0 **>+l 

Expanding this equation in h, we obtain for h > 1 



ij 2 



'd.2+1 



2Q(x(j>)) (4-10) 



J4-11) 

We consider this equation when p is near a branch point multiply it by | J^f^) > 
take the residues when p —>■ fi a and sum over all the branch points. Notice that 
| y^^y is finite at p = p, i.e. it has no poles at branch points. 

Let us first compute the LHS: 



V Res 

Z — ' P~>Mo 

V Res 

V Res 

a 

V Res 



y(p) - y(p) 

dE Pif (q) \y(p)W%(p) + y(p)Wj$(p) + E^vCp^JV) 



y(p) - y(p) 

dE p>¥ (q) y(p)W^(p)+y(p)W^(p) 



yip) - yip) 

dE prp (q) \v(p)W$$(jp) - y{p)(W^(p) + Y^W^ltf)) 



yip) - yip) 



V Res 



!/(p)<(p) -y(pW,o(p) 



yip) -yip) 

^Res d^(g)Wi ( J?(p) 

* ' p— >/Li Q 

]T Res dS p , (g) W$(p) - d%, (g) W$(p) 

* — ' p— >/i a 
a 

Res ^ Pi0 (g) ^'(p) - W[%)) 
£ ^s d5 Pi0 (g) (2Wi<J>(p) + £ W$(p<)) 

p l ^p,p 

2^ Res ^(g)^»(p) 

r- E (7 ^fJW / fi (p>?) + / <W / 5 G^)) 

^ \^v4i ^-4i / 



-2 Res ^(g)^^) 
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= 2Wf } (g) 
(4-12) 

Let us compute the RHS in a similar way. It reads 



V Res 



{yip) - y(p))dx(p) 



i=0 



/i-l 



+ E <o VK'T'V) + 2QW(x(p))dx(p) 



m=l 



V Res 

J—^ r>->-u„ 



P^na {y{p) - y{p))dx(p) I 



w% 1 \p,p) + E <?(p)<o~ m) (p) 



+w^- 1) (p,p) + e w$<p)w§r*<p) 



h-l 



m=l 



2 V Res 

f ^ r) — ► !/.„ 



\dE p ip(q) 



wf; 1] { P , P ) + E <o } (p)<- m) b) 



m,=l 



/i-l 



m=l 



>->Ma {y{p) - y{p))dx(p) 

(4 - 13) 
Thus we have, for /i > 1: 

<» m = y Res ij^MKTW) + ^ <"w<"'"(rt) (4 . 14) 



(y(p) - y{p))dx{p) 



Using properties (|4-5jl and (|4-7jl . this last equation can be changed into (we add 
terms which don't have poles at branch points, and thus whose resiudes vanish): 

wi% } = - y Res ^)(<- 1 Vrt + s^,^'(P)^r ) w) (4 . 15) 



(y(p) - y{p))dx(p) 



This is the case = of the following theorem (proved below by recursively applying 
the loop insertion operator d/dVi): 

Theorem: 



~ E Q Res 



P-*Mq (y(p)-y(p)) dx(p) 

( m ) f„ ^ ^ u/( /i_m ) ^ 



(4-16) 



which can be diagrammatically represented by 





(4-17) 



where a sphere with h holes and k legs is the fc-point function to order h wj: h \ and 
the arrow means the following integration: 

/ q -\dE q Jp) 

P — ► K = Y Res rrr /4 j ? n 4 -i8 

This equation is a recursion relation, whose initial condition is given by (|2-33|) . i.e. 

B(p,q). (4-19) 
proof: 

We have already proved theorem ()4-16|) for k = 0. Let us now show that the rule 
()4-18|) is compatible with the derivation wrt the potential V\. 

One preliminary needed formula is the action of the loop insertion operator on the 
function y(p). It is well known (see for exemple ^H]) that: 

dM dx(r) = -^4 (4-20) 




dVi(x(r)) dx(p) 
First, the action of the loop insertion operator on the Bergmann kernel gives 12 II : 

ggfog) dx(r) = V Res dE a^ B (P~0 B ^ r ) 

Reg 1 dE a (q) [Bjp, flflfc r) + B(r, QB(£, pj\ 

(4-21) 

Note that this corresponds exactly to the conjecture for k = 1 and h = 0. 

By integrating this expression with respect to p, we obtain the action on the arrowed 
edges: 

ddE prp (q) = v Reg 1 dE a (q) [dE prp (QB(£,r) + B(r,QdE prp (0] 

WMr)) V^"- 2 (y(O-y(O)^(O 

where q, p and r are outside the integration contour. 

There is one last quantity to evaluate, corresponding to the vertices: 

9 ( 1 7^) dx(r) = , B ^ r) -fSf) , (4-23) 



dVi(r) \y(p) - y(p) J (y(p) - y(p)) 2 dx(p) 

Let us now interpret diagrammatically these relations. 

Eq. (|4-21j) represents the action of the loop insertion operator on the non arrowed 
edges: 




d 



p 

(4-24) 
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In order to interpret ([4-22)1 . one has to move the point p inside the integration 
contour for £. For this purpose, let us consider the action of the loop operator on the 
vertex, dx(r) av * x{r)) £ Q Res^ dE % qip) m) _ y \ q))dx{q) • It gives the contribution 



V Rps dEg^jp) B(q,r)-B(q,r) 

Z-^a rusB 2 (y(q)-v(a)) 2 dx(o 



{y(q)-y(<i)) 2 dx(q) 2 _ oc .n 

^(p) dE g , 7 (QB($,r)+rfg g , 7 (g)B(g,r) l 4 "^ 



where q lies outside the integration contour for £. We move the integration contour for 
£ in the second term so that g lies inside. One has poles only when £ — > q and £ — > g 
which can be written 

Res Res = Res Res — Res Res — Res Res . (4-26) 

q^Va t;—>fj.ct f— >iJ. a q-*^a q^t^a £,^q q^y-a ^q 

These contributions cancel totally the first term in ([4-25)1 . Thus one finally obtains: 

dx (r) - V Res dEqM - = 

1 ) dV^xir)) 2 (y(q) - y{q))dx{q) 

= V Res Res - _ dE *&) dE qrq (QB(£, r) + dE qrq (QB(l r) 
^ 2(y(£) - y{S))dx(() 2{y{q) - y{q))dx{q) 

(4 - 27) 

which corresponds to adding one leg to the vertex: 



dx(r) smm > ^ 



+ ) 4 ) < 



dVi(x(r)) ' ^ " — )— 

(4-28) 

This ensures that the theorem is true. □ 

Therefore we have the following diagrammatic rules: 

Theorem 4.1 The k + 1 point function to order h, is the sum over all possible 

diagrams obtained as follows: 

• choose one of the variables, say p\ as the root. 

• draw all rooted skeleton binary trees (i.e. trees with vertices of valence < 3) with 
k — 1 + 2h edges. Draw arrows going from root to leaves, that puts a partial 
ordering on vertices. A vertex V\ preceeds V 2 if there is an oriented path from V\ 
to V2, and a vertex always preceeds itself 

• add, in all possible ways, k non arrowed edges ending at the point p%. . . ,Pk+i, 
so that each vertex has valence < 3. 

• complete the diagrams in all possible ways with h non arrowed edges joining 2 
comparable (with the partial ordering of the tree) vertices, so as to get diagrams 
with valence 3 only. 
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• At each vertex, mark one leg with a dot, in all possible inequivalent ways. In 
general there are 2 possibilities at each vertex, except if the 2 branches coming 
out of the vertex are symetric (either 1 non-oriented edge, or 2 identical subgraphs 
with no external legs). Each inequivalent diagram is counted exactly once (this is 
the standard way of computing symmetry factors). 

• assign to each such diagram a value obtained as follows: each non-arrowed edge 
is a Bergmann kernel, a residue of the form \4-l£ty is computed at each vertex, 
where the q variable corresponds to the marked edge. The ordering for computing 
residues is given by the arrows, from leaves to root. 

Notice that those rules are extremely similar to those first found for the 1-matrix 
model in The 1-MM is merely the reduction of those rules to the case of 

an hyperelliptical surfaces. 

5 Topological expansion of the free energy 

In order to find the topological expansion of the free energy, one has to "integrate 
wrt the potential V\ n the expansion of the one point function. First let us define the 
H x and H y operators corresponding respectively to the "inverse" of the loop insertion 
operator ^ and 

5.1 The H operators 

We define the operators H x and H y such that, for any meromorphic differential form 
0, one has 

H x .(f> := Res V^x) - Res (V 2 (y) -xy)<j) + T [ " + Y] e f / (5-1) 

°°* °°y Joox : Jb, 



and 



roo x r 

H y 4 := Res V 2 (y) - Res (Vi(x) — xy) (f> + T <p - V e { <b 0. (5-2) 

JoOy V ^ 

Note that 

(H x + H y ).(p= Res xy.cj) (5-3) 
Note that if = df is an exact differential, we have: 

Hx-df = — Res f ydx , H y .df = — Res f xdy. (5-4) 

We now compute the action of H x on the two-point function on the Bergmann 
kernel B(p, q) (in this computation H x acts on the variable p), that is 

H x .B(p, q) = Res Vi(x(p)) B(p, q) - Res (V 2 (y(p)) - xy) B(p, q) 
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+T / V B{p,q) + ^i <f B(p,q). 

(5-5) 



We integrate by parts the residues at infinities. Provided that B{p,q) = d p dS p , (q), 
that V{(x(p)) - y{p) ~ as p -> oo x . and that V£(y(p)) - x(p) ~ as p -> oo„, we 



can write 

T 

H x .B(p,q) = - Res — -dx(p) dS P:0 (q) - Res y(p)dx(p) dS Pt0 (q) 

oo x X{p) co x ,co y 

+ Res —-^-dy(p) dS PtQ (q) + TdE^ yt00x {q) + Y] $ B(p, <?). 

y(p) i M 

(5-6) 

Finally, the Riemann bilinear identities [2131211] associated to the fixed filling fractions 
conditions give 

H x .B(p, q) = TdSoa (q) - Res y(p)dx(p) dS Pt0 (q) - Res y{p)dx(p) dS Pj0 (q) 

OQ X OOy 

-T dS^q) + TdS^M + I> ^ B ^ 
= Res y{p)dx dS p a (q) 

1 



+ Y j e i l B(p,q) 

= -y{q)dx{q). (5-7) 

That is: 

H x .B(.,q) = -y(q)dx(q). (5-8) 

Remark 5.1 The leading order of the free energy is already known 1 : 

2 jM = Res (Vl (x) + V 2 (y) - xy) ydx + T f ydx + Y e { <L ydx 

J OOy V JBi 



(5-9) 
In other terms, one has 



2T {0) = H x {ydx) (5-10) 



5.2 Finding the free energy 

Theorem 5.1 For any h we have : 

{\-2h)W§=H m .W%$ (5-11) 
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Notice that this is already proven for h = 0. 
proof: 

Let us study how H x acts on the diagrams composing W^j ■ By symmetry one can 
consider that H x acts on the leaf of the diagram 2 . Two different configurations can 
occur. 

If the leaf is linked to an innermost vertex, one has to compute 

dE r>T (q)H x .B(r,.)B(r,p)f(p) ^ dE r , r (q)y(r)B(r,p)f(p) 



E Reg a^r,ryg)^x-J=>{i,-)^{i,p)j \Pj = _ Reg 



Jr^n* 2((y(r) - y{r))dx{r) ^ 2 ((y( r ) ~ v( r )) 

0. 

(5 - 12) 

Otherwise, H x acts on B(r, .) where there is an arrow going out of the vertex 
corresponding to the integration of r: 

A = y Res V Res dE r ^{q)dE p ^r)y{r)f{p) 

^ r^ a , P ^ a 2(y(r) -y{r)){y{p) - y{p))dx(r)dx(p) 

By moving the integration contours so that r lies inside the contour of p, one keeps 
only the contributions when r — > p and r — > p, that is 

A = y Res R es dE rt r(q)dE p>v (r)y(r)f(p) 

a*"** r ^ p 2 (^( r ) - y( f ))(y(p) ~ y{p))dx{r)dx{p) 
+ V Res Res dE r,r(Q) dE P,p( r )y( ¥ )fip) 

a^ 1 ""* 2 (y( r ) -y( f ))(y(p) - y{p))dx(r)dx{p) 

= r Res (5-14) 
^ 2(y(p) - y(p))dx(p) 

The same diagram with H x acting on B(r, .) gives the same contribution and cancels 
the | factor. 

Diagrammatically, this reads 



+ = i i (5-15) 



and 



H x 



0. (5-16) 



2 The diagrams composing q have two external legs, one root with an arrow and one leaf which 
is a Bergmann kernel. 
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Thus, because there are 2h—l arrowed edges composing any graph contributing to 
one has: 

= -(2h-l)wi$(p). (5-17) 

□ 

Beside, we have the following property: 
Lemma 5.1 For any h^l: 

H x .wW(q) - H y .W^(q) = (5-18) 
whose rather technical proof is found in appendix B. 

Corollary 5.1 One easily derive from this theorem that for any h ^ 1, the free energy 
can be written as 



(2 - 2h)jM = -H X W% = -H y W$, 



(5-19) 



up to an integration constant which does not depend on V\ or V 2 . 
5.3 Dependence on other parameters 

Now, with the explicit dependence of the free energy on the potentials V\ and V 2 
in hands, one is interested in the dependence in the other momenta, i.e. the filling 
fractions and the temperature T. 

The definitions of the filling fractions implies that 



dy{p)dx{p) 



= 2mdui= i B(p,q), (5-20) 

JBi 



where the du^s denote the normalized holomorphic differential on cycles A4. This 
allows to check that the derivation wrt the filling fractions is compatible with our 
diagrammatic rules through 

dB(p, q) jf d 
dei 

which generalizes to 



lidm BiM) - (5 - 21) 

^ = /<W (5-22) 



In the same fashion, one derives 

or 



/ W$>(p,q). (5-23) 
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5.4 Scaling equation 

For the sake of completeness and in order to interpret the formula found for the free en- 
ergy, one has to consider a slightly more general model, in which the coupling constant 
between the two matrices is not fixed, i.e. we consider the partition function 

Z:= I dMldM2e -\<VrmM.)+V 2 {Mz)-nM,M,) = -T ^ (5.34) 

where k is an arbitrary non vanishing coupling constant. 

One can easily obtain this model by rescaling our momenta: 

tk — > — , tk — > — and T — > — . (5-25) 

Kj Kj K 

After this redefinition of the parameters, all the results presented in this paper 
remain valid 3 and one obtains 

(2 - 2h)J^ = -H x w[fi + f(T, <*, k), (5-26) 

with 

H x .<j> := Res Vi(x) <\> - Res (V 2 (y) - Kxy) <p + T [ " + V / (j). (5-27) 

Joo x V hi 

Now, define 

\/h ^ 1 FW = ; . (5-28) 

Eq. (jEZZD and ((E2BD give 

dFW dH x W^(p) 



(2 - 2h)- 



dei dei 



WW{p)-H x - 



(5 - 29) 

That is to say 



-(2-2/,)/ W%(p). 



^ = - / W[%) and ^ = " / ' < } (P)- (5-30) 



On the other hand, by definition 

dt k ~T ' lv " Vjryy (9Vi(a;(p)) 



fc=l 



3 Note that the diagrammatic rules remain the same, only the Riemann surface is affected by this 
rescaling. 
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Which leads to 

dl QpW 



k=l 

One can also derive 



I>^^ = - Res VMp^W^ip). (5-32) 



2>^r = -% s 



it=i 



oo„ (2-2/i) dF 2 (y(p)) 



4 s (2-2/1) av 2 (y(p)) 

+ 4 S (2-2/i) 9y 2 (j/(p)) 
Res ^(p))!*^). 



(5 - 33) 



The last identity uses the property (|5-3|) of the H operators. 

There is one additive momentum that one has to consider now, which gives 

Qp(h) 

— — = - Res x(p)y(p)W 0A (p) (5-34) 

OK OOy 

Thus, putting everything together, for any set of filling fractions e$, any coupling 
constant k and any temperature T, fulfills the homogenous equation 



A? — 1 ft! — 1 2 



(5-35) 

Hence, is our final answer for the free energy: 

p{h) = jr{h) = _J_ H W W = _J_ H w (h) (5 _ 36) 
2/i-2 W 2h-2 y C ' 1 ) v ; 

Remark 5.2 One can note that l)5-35j) is nothing but saying that rescaling the expansion 
parameter H is equivalent to rescaling all the other momenta in the same way: 

8F ^ 8F %^~dP 8F dF dF , „ N 

° = % + + E ^ + ^- + k- + £ (5-37) 

fc=i fc=i 2 

5.5 Explicit computation of the free energy: diagrammatic 
rules 

The equation ()5-36|) allows to compute any term in the expansion of the free energy, but 
the first correction to the leading order, using diagrammatic rules. For this purpose, 
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one just introduce a new bi-valent vertex corresponding to the action of the H x operator 
on the root of the diagrams composing the one point functions. One represent it, 



/ Res / / \ — , 7 x (5-38) 
^™ {y{q) - y{q))dx{q) 



Using this new vertex, one obtains the Tt'th order correction to the free energy by 
considering the sum of the diagrams contributing to Wiq where the ingoing arrow is 
replaced by this new bi-valent vertex. 

Remark 5.3 We must remind that this technic does not give the first order correction. 
Nevertheless this particular correction has already been computed in the literature |2Uj . 



6 Conclusion 

In this paper, we have obtained a closed expression for the complete expansion of the 
free energy of the formal hermitian two-matrix model as residues on the spectral curve 
giving an answer to one of the problems left undone in p2 1 . On the way, we also 
refined the diagrammatic rules used to compute non-mixed correlation functions. We 
particularly exhibit that these functions actually only depend on what happens near 
the branch points. The link with the diagrammatic technique of JH] for the 1MM looks 
also more evident. 

Nevertheless, there still remain many different problems one should address from 
this starting point. In this paper, we begin the computation of correlation functions 
mixing traces in M\ and traces in M2, but the general case does not seem as simple 
for the moment. In the same way, using the result of [17] . one would like to compute 
the expansion of mixed traces. 

This natural generalization of diagrammatic techniques from the 1MM to the 2MM 
also points out that the study of the chain of matrices |2Ij would give a more general 
knowledge of the link between algebraic geometry and matrix models. 

Let us also mention that we have considered here only polynomial potentials, but 
it is clear that the whole method should extend easily to the more general class of 
semi-classical potentials 011], whose loop equations are very similar 23 j. The 1-MM 
with hard edges was already treated in [T^] . 

Another extension of this problem is to compute the 1/N expansion of random 
matrix models for non-hermitian matrices. This problem is already partialy addressed 
in [SHHHUIHE] , but a diagrammatic technique has not been found yet. 

Finally, as we already mentioned it, solving the formal model is different from 
solving the physical problem corresponding to a convergent integral and constrained 
filling fractions lowering the energy. This problem, which is expected to be similar 
to jH] is not addressed for the moment and should give a natural complement to this 
work. 
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Appendix A Computation of Wj^. 

In this section, we compute a new set of correlation functions including one trace in 
M 2 . In other words, we compute all the wj* s. More precisely, we explain how the 
diagrammatic rules are modified when one changes one trace 

Tr Tr . (1-1) 

x(p) — Mi y(p)-M 2 

We do not give explicit formulas for all the with I > 1. This problem will be 
addressed in another work. 

A.l Examples 

In order to give some idea of the result, let us review some already known functions [T9l 

Emu. 

A. 1.1 2-point functions 

W$(p, q) = W${p, q) = -W±°}(p, q) = B(p, q) (1-2) 

A. 1.2 3-point functions 

m/(°)/ \ t? B (puq) B (P2,q)B(p 3 ,q) ( q\ 

w »<*^=?jz — mm — (1 " 3) 

W (0) (p 1 p 2 Pz) = - V s Res ^(Pi'g)^^,?)^^,?) _ Reg B(p 1 ,q)B(p2,q)B(p 3 ,q) 
" l ' ' ' g-^Mc dx(q)dy(q) g-»P3 dx(q)dy(q) 

(1-4) 

proof: 

d 

W 2 ,i(pi,P2,q) = a , r( - ( A w ljl (p l ,q))dy(q)dx(p l )dx(p 2 ) 

dVi{x(p 2 )) 

d dxicf) 
~ x\r( I ^ ( w 2,o(Pi, q)-f T -My{q)dx{p 1 )dx{p 2 ) 
dVi(x(p 2 )) dy(q) 

W f \ I W 2 fl{px,q), fW 1A (p 2 ,q) 

= -W 3fi {pi, p2, q) + — j j~t — a,' 



dx(q) V d y(q) 
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| W 1 , 1 (p 2 ,q) d f W 2 ,o(pi,q) \ 
dy(q) q V dx(q) J 

= -Ws,o(P^ Q) + d q { d^iq) ) 

(1-5) 

□ 



A. 2 General case 
Theorem A.l For any \K\ + h > 1 

W$l 1 (pi,---,Pk,q)+W$ l+lj0 (p l ,...,p k ,q) = d q f ih \pi,...,p k ,q) (1-6) 

where f^ h \px, ■ ■ ■ ,Pk, q) is obtained from the diagrams composing (pi, . . . , P k) by 
cutting any of its propagator as follows: 

P . P 2-Z^jW P l q q P 2 ^ 

and 

- mw> p q q ^~ ( ' (1 " 8) 

proof: 

W$ ^Pk, q) can be obtained by differentiating W$ (pk) wr t V2(y(q)). So one has 
to know the action of on the diagrammatic rules. 
One can reexpress the 3-point function of (jl-4j) as 

9 Dt v _ 9 ( B(s,p)B(s,i 

dVM*)) [<1,P) ~ MM*)) { } V dx(s)dy(s) 
This can be graphically represented by 




d d 

p i = - nr, i \ \ p q 



dV 2 (y(s)) r 9V 1 (x(a)) 

s \dx(s)dy(s) ^ S 
(1-10) 

On the other hand, by integration, this implies that 

9 , - _ dE a: g(q) [dE Pir (a)B(a, s) + B(s,a)dE p>r (a)] 

dV 2 (y(s)) pAq) 2(y(a) - y(a))dx(a) 

- Res Res dE aA<l) dE pA b ) B ( a i s) 

^ a-^na b^a (y(b) — y(a))(x(b) — x(a)) 

(1-11) 
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Thus the derivative wrt V 2 of the first vertex, » T/ 9 , Res„^,, „, ,-, Ep, f^), , , 

' oV 2 {y{s)) ^~>a P^M-a 2(y(p)-y(p))dx(p) 

gives: 

dK,a(?)[d£pjK") B ( a > s )+ B ( s > s ) dE P,p( a )] 
" Z^a 2^/3 ^ es p^Ma -K-eSa^M/J 4(y(a)-y(a))dxfaXy(pVj/(p))d:r(p) 

-V Res Res Res,. dE a Jg)dE p , ¥ (b)B(a,s) (1-12) 

n,tJb p^ Ma ixeb a _» s ru?b 6^ 2(y(b)-y(a))(x(&)-x(a))({/(p)-j/(p))dx(p) ' V y 

_i_ V Rps dE P,p(<l) d 1 

^Z^a rLtJb P- > Ma 2 aV 2 (s) (j/(p)-»(p))rfx(p) 

where p lies outside the integration contour for a. Note that all the derivations have 
been performed by keeping the x's fixed. 

Moving the integration contours so that p lies inside gives a contribution that cancels 
the last terms and one obtains 



V Res dE P,p(i) 8 v Res 

dV 2 (y(s)) 2-~<ol J^^p-^/ia 2(y(p)-y(p))dx{p) <9Vi(x(s)) lxc& P-^a 2(j/(p)-s/(p))cia;(p) 



dE p> -p(q) 



" s ^ dx(s)dy(s) ^^p^^ a 2(y(p)-y{p))dx(p) J 

(1-13) 



This graphically reads: 
d , / d 



dV 2 (y(s)) n dViMs)) 

1 



-cL 



dx(s)dy(s) ^ 



(1 - 14) 

These two relations prove straightforwardly the theorem. □ 

Remark A.l For \K\ = 1, one can go further in the description of (p±,p 2 ) + 

W^Q (pi,p 2 )- Indeed, the graphs obtained by cutting the propagators of W^ipi) have no 
pole when p\ — > ^2 except when one cuts the root of the graph, in which case it gives the 
contribution ^g^) W-ffffo), that is to say: 

wffcpi.K) + = -d P2 {^^y<ofe) + /(pi >P2 )} , (1-15) 

where f(pi,p 2 ) has no pole as pi —> P2- 

Appendix B Symmetry of the free energy under the 
exchange of x and y 

For h ^ 1 

" = {H Xp H yp2 - H yp H Xpi ) W^(p u p 2 ) (2-1) 

The only terms that contributes in the exchange of H x and H y come from the 
integration contours around oo x and oo y , i.e. 

A {H Xpi H yp2 - H yp2 H Xpi ) W[ h l ( Pl , p 2 ) 



27 



= - Res Res V l (x(pi))(V l (x(p 2 )) - x(p 2 )y(p 2 ))Wf{(p 1 ,p 2 )- 
- Res Res V 2 (y(p 2 ))(V 2 (y(pi)) - x(pi)y(p 1 ))W^(p 1} p 2 ). 

(2-2) 



W[% u p 2 ) = -WgfCpi.pa) - rf P2 (t^It^SW + /(pi,J*) } • (2-3) 



On the other hand, one knows from remark (jl-15j) . that 

[ B(p 1 ,p 2 ) 
\dx(p 2 )dy(p: 

where f{pi,p 2 ) has no pole when p 2 — > p\. Note that W 2 q (pi,P2) do not have pole 
when p 2 — > pi either. 

Thus, the only non vanishing terms in (|2-2|) come from d P2 

This reads 

Res Pl ^ x Res P2 _ Pl V 1 (x(p 1 ))(V 1 (x(p 2 )) - x(p 2 )y(p 2 ))d P2 {^yWjJW} 

+ Hes^oo, Res P2 _ Pl V 2 (y(p 2 ))(V 2 (y( Pl )) - x(p 1 )y(p 1 ))d P2 { j^% 2) W^{p 2 ) 

(2-4) 

We now evaluate these residues by part and obtain that 



A = - Res V 1 (x(p))d p 

p^co x 



(V{(x(p)) - y{p))dx(p) - x{p)dy(p) ' 



- Res (V 2 (y(p)) - x(p)y(p))d p 

JJ^OOy 



dx(p)dy(p) 

vMp)W^(p) 

dx{p) 



(2-5) 

Another integration by part can be written 

A = + Res ^(pjj ^^^-^Pfy-^^V wCp) 
p-^oo x dy(p) 

V 2 (y(p))W^(p) 

+ Res {V 2 (y(p))dy(p) - x(p)dy(p) - y{p)dx(p)) ' . 

P^OOy dX{P) 

(2-6) 

Knowing that V{(x(p)) - y{p) ~ ^ when p ^ oo x and V 2 \y(p)) - x(p) ~ ^ 
when p — > oo. y and that the other factors do not have poles at infinities, one can finally 
write that 

(H Xn H yp2 -H yp2 H Xpi )W%( Pl ,p 2 ) = - Res x{q)y{q)W[%). (2-7) 

Let us now show that the RHS vanishes by looking precisely at the form of }{q)- 
Because it comes from a vertex, ()4-16j) . with moving the integration contour, implies 
that 

2 Res x(q)y(q)W^(q) = -2V Res x{q)y{q)W[%) 



q— >^a 
a 
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V V Res Res ^M^pM x 
p^ M/3 _ y(p))dx(p) 



a a 



X 



(2-8) 



h-l 



w& _1 W) + E <V)< m) (p) 



m=l 



Now we move one more time the integration contour and clearly segregate the case 
when p and q are near the same branch point. 



2 Res x{q)y{q)w[%) 



a^beta 



x(q)y(q)dE P: p(q) 

""'^ -y(p))dx(p) 



ERes Res 
T) ► //.„. fl -tll.r. 



X 



h-l 



+ E <(ri< m) (p) 

m=l 

x(q)y(q)d~E P: p(q) 



Res Res — Res Res_ 

p^fi a q->fi a q^fJ. a p^q , q 



X 



(2-9) 



{yip) - y(p))dx(p) 
w^W) + E <b)< m) b) 



h-l 



m=l 



One can see that the integrand has no pole as q approaches a branch point. Thus 
one keeps only 



ERes Res_ 
a — w/.„, r> — \a n 



x(q)y(q)dE p! p(q) 



h-i 



'</-/*« p-9. a (y(p) -y{p))dx(p) 
We finally perform the integration and get 

2 Res x{q)y{q)W[%) = E Res TTT 
00^,00,, ' q^n a dx(q) 

= 0° 



K~ 1] (p,p) + E <W< m) b) 



m=l 



• (2-10) 



h-i 



< _1 W) + E <(f)< m) w 

(2-1 



m=l 



QED. 



References 

[1] G.Akemann, "Higher genus correlators for the Hermitian matrix model with mul- 
tiple cuts", Nucl. Phys. B482 (1996) 403, hep-th/9606004 

[2] G.Akemann and J.Ambj0rn, "New universal spectral correlators", J. Phys. A29 
(1996) L555-L560, cond-mat/9606129. 

[3] J.Ambj0rn, L.Chekhov, C.F.Kristjansen and Yu.Makeenko, "Matrix model calcu- 
lations beyond the spherical limit", Nucl. Phys. B404 (1993) 127-172; Erratum 
ibid. B449 (1995) 681, hep-th/9302014. 



29 



[4] M. Bertola, "Free Energy of the Two-Matrix Model/dToda Tau-Function", 
preprint CRM-2921 (2003), hep-th/0306184. 

[5] M. Bertola, "Second and Third Order Observables of the Two-Matrix Model", J. 
High Energy Phys. JHEP11 (2003)062 , hep-th/0309192. 

[6] M. Bertola, "Bilinear semi-classical moment functionals and their integral repre- 
sentation", math/02051. 

[7] M. Bertola, "Two- matrix model with semiclassical potentials and extended 
Whitham hierarchy", hep-th/0511295. 

[8] G. Bonnet, F. David, B. Eynard, "Breakdown of universality in multi-cut matrix 
models", J. Phys. A33 6739-6768 (2000). 

[9] E. Brezin, C. Itzykson, G. Parisi, and J. Zuber, Comm. Math. Phys. 59, 35 (1978). 

[10] L.Chekhov, "Genus one corrections to multi-cut matrix model solutions", Theor. 
Math. Phys. 141 (2004) 1640-1653, hep-th/0401089. 

[11] L.Chekhov and B. Eynard, "Hermitian matrix model free energy: Feynman graph 
technique for all genera", hep-th/0504116 

[12] L.Chekhov, "Matrix models with hard walls: Geometry and solutions"; contribu- 
tion to special volume of J. Phys. A on matrix models, hep-th/0602013. 

[13] David F., "Loop equations and nonperturbative effects in two-dimensional quan- 
tum gravity". Mod. Phys. Lett. A5 (1990) 1019. 

[14] P. Di Francesco, P. Ginsparg, J. Zinn- Justin, "2D Gravity and Random Matrices", 
Phys. Rep. 254, 1 (1995). 

[15] R.Dijkgraaf and C.Vafa, "Matrix Models, Topological Strings, and Supersymmet- 
ric Gauge Theories", Nucl.Phys. 644 (2002) 3-20, hep-th/0206255; "On Geometry 
and Matrix Models", Nucl.Phys. 644 (2002) 21-39, hep-th/0207106; "A Pertur- 
bative Window into Non-Perturbative Physics", hep-th/0208048. 

[16] B. Eynard, "Topological expansion for the 1-hermitian matrix model correlation 
functions", JHEP/024A/0904, xxx, hep-th/0407261. 

[17] B. Eynard, N. Orantin, "Mixed correlation functions in the 2-matrix model, and 
the Bethe ansatz", J. High Energy Phys. JHEP08(2005)028, hep-th/0504029. 

[18] B. Eynard, "Large N expansion of the 2-matrix model", JHEP 01 (2003) 051, 
hep-th/0210047. 

[19] B. Eynard, "Large N expansion of the 2-matrix model, multicut case", preprint 
SPHT03/106, ccsd-00000521, math-ph/0307052. 

[20] B. Eynard, A. Kokotov, and D. Korotkin, U l/N 2 corrections to free energy in 
Hermitian two-matrix model", hep-th/0401166. 



30 



[21] B. Eynard, "Master loop equations, free energy and correlations for the chain 
of matrices", J. High Energy Phys. JHEP11(2003)018, hep-th/0309036, ccsd- 
00000572. 

[22] B. Eynard, "Polynomes biorthogonaux, probleme de Riemann-Hilbert et geometrie 
algebrique", Habilitation a diriger les recherches, universite Paris VII, (2005). 

[23] B. Eynard, "Loop equations for the semiclassical 2-matrix model with hard edges", 
math-ph/0504002. 

[24] B. Eynard, N.Orantin, "Topological expansion of the 2-matrix model correlation 
functions: diagrammatic rules for a residue formula", J. High Energy Phys. 
JHEP12(2005)034, math-ph/0504058. 

[25] H.M. Farkas, I. Kra, "Riemann surfaces" 2nd edition, Springer Verlag, 1992. 

[26] J.D. Fay, "Theta functions on Riemann surfaces", Springer Verlag, 1973. 

[27] V.A. Kazakov, "Ising model on a dynamical planar random lattice: exact solu- 
tion", Phys Lett. A119, 140-144 (1986). 

[28] V.A. Kazakov, A. Marshakov, "Complex Curve of the Two Matrix Model and its 
Tau-function", J. Phys. A36 (2003) 3107-3136, hep-th/0211236. 

[29] I.K.Kostov, " Conformal field theory techniques in random matrix models", hep- 
th/9907060. 

[30] I.Krichever "The r-function of the universal Whitham hierarchy, matrix models 
and topological field theories", Commun.Pure Appl.Math. 47 (1992) 437; hep- 
th/9205110 

[31] L.Chekhov and A.Mironov, "Matrix models vs. Seiberg-Witten/Whitham theo- 
ries", Phys.Lett. 552B (2003) 293-302, hep-th/0209085 

[32] L.Chekhov, A.Marshakov, A.Mironov, and D.Vasiliev, "DV and WDVV", Phys. 
Lett. 562B (2003) 323-338, hep-th/0301071 

[33] M.L. Mehta, Random Matrices ,2nd edition, (Academic Press, New York, 1991). 

[34] M. Staudacher, " Combinatorial solution of the 2-matrix model" , Phys. Lett. B305 
(1993) 332-338. 

[35] G. 't Hooft, Nuc. Phys. B72, 461 (1974). 

[36] P. Wiegmann, A. Zabrodin, 'Large N expansion for normal and complex matrix 
ensembles', hep-th/0309253. 

[37] I. Krichever, M. Mineev-Weinstein, P. Wiegmann, A. Zabrodin, 'Laplacian Growth 
and Whitham Equations of Soliton Theory', nlin. SI/0311005. 

[38] A. Zabrodin, P. Wiegmann, 'Large N expansion for the 2D Dyson gas', hep- 
th/0601009. 



31 



